Temperature dependence of thermal conductivity in 1D nonlinear lattices

We examine the temperature dependence of thermal conductivity of one dimensional nonlinear (anharmonic) lattices with and without on-site potential. It is found from computer simulation that the heat conductivity depends on temperature via the strength of nonlinearity. Based on this correlation, we make a conjecture in the effective phonon theory that the mean-free-path of the effective phonon is inversely proportional to the strength of nonlinearity. We demonstrate analytically and numerically that the temperature behavior of the heat conductivity $\kappa\propto1/T$ is not universal for 1D harmonic lattices with a small nonlinear perturbation. The computer simulations of temperature dependence of heat conductivity in general 1D nonlinear lattices are in good agreements with our theoretic predictions. Possible experimental test is discussed.Comment: 6 pages and 2 figures. Accepted for publication in Europhys. Let

Global Birkhoff coordinates for the periodic Toda lattice

In this paper we prove that the periodic Toda lattice admits globally defined Birkhoff coordinates.Comment: 32 page

Bifurcations of discrete breathers in a diatomic Fermi-Pasta-Ulam chain

Discrete breathers are time-periodic, spatially localized solutions of the equations of motion for a system of classical degrees of freedom interacting on a lattice. Such solutions are investigated for a diatomic Fermi-Pasta-Ulam chain, i. e., a chain of alternate heavy and light masses coupled by anharmonic forces. For hard interaction potentials, discrete breathers in this model are known to exist either as ``optic breathers'' with frequencies above the optic band, or as ``acoustic breathers'' with frequencies in the gap between the acoustic and the optic band. In this paper, bifurcations between different types of discrete breathers are found numerically, with the mass ratio m and the breather frequency omega as bifurcation parameters. We identify a period tripling bifurcation around optic breathers, which leads to new breather solutions with frequencies in the gap, and a second local bifurcation around acoustic breathers. These results provide new breather solutions of the FPU system which interpolate between the classical acoustic and optic modes. The two bifurcation lines originate from a particular ``corner'' in parameter space (omega,m). As parameters lie near this corner, we prove by means of a center manifold reduction that small amplitude solutions can be described by a four-dimensional reversible map. This allows us to derive formally a continuum limit differential equation which characterizes at leading order the numerically observed bifurcations.Comment: 30 pages, 10 figure

Riemann solvers and undercompressive shocks of convex FPU chains

We consider FPU-type atomic chains with general convex potentials. The naive continuum limit in the hyperbolic space-time scaling is the p-system of mass and momentum conservation. We systematically compare Riemann solutions to the p-system with numerical solutions to discrete Riemann problems in FPU chains, and argue that the latter can be described by modified p-system Riemann solvers. We allow the flux to have a turning point, and observe a third type of elementary wave (conservative shocks) in the atomistic simulations. These waves are heteroclinic travelling waves and correspond to non-classical, undercompressive shocks of the p-system. We analyse such shocks for fluxes with one or more turning points. Depending on the convexity properties of the flux we propose FPU-Riemann solvers. Our numerical simulations confirm that Lax-shocks are replaced by so called dispersive shocks. For convex-concave flux we provide numerical evidence that convex FPU chains follow the p-system in generating conservative shocks that are supersonic. For concave-convex flux, however, the conservative shocks of the p-system are subsonic and do not appear in FPU-Riemann solutions

Continuous Symmetries of Difference Equations

Lie group theory was originally created more than 100 years ago as a tool for solving ordinary and partial differential equations. In this article we review the results of a much more recent program: the use of Lie groups to study difference equations. We show that the mismatch between continuous symmetries and discrete equations can be resolved in at least two manners. One is to use generalized symmetries acting on solutions of difference equations, but leaving the lattice invariant. The other is to restrict to point symmetries, but to allow them to also transform the lattice.Comment: Review articl

Heat flux distribution and rectification of complex networks

10.1088/1367-2630/12/2/023016New Journal of Physics12